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Abstract: Heavy quark production in 77 collisions is reanalyzed. It is argued that evaluating 
the cross section (7(77 — ► QQ) in a well-defined renormalization scheme requires the inclusion of 
direct photon contributions up to the order a 2 a 2 . The order a 2 a 2 , direct photon contributions 
are furthermore needed for factorization scale invariance of the sum of direct and resolved photon 
contributions. The importance of quantitative analysis of renormalization and factorization scale 
dependence of the approximation currently used for the evaluation of 17(77 — * QQ) is emphasized 
as the only way of estimating the theoretical uncertainty related to the ambiguity in choosing 
these scales. 
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1. Introduction 



Heavy quark production in 77 collisions has recently received increased theoretical attention 
[1, 2, 3] motivated by new experimental data on cc and bb production from LEP2 [4], which 
provide particularly suitable framework for the confrontation of perturbative QCD with data. 
The results obtained so far and shown, for instance, in Fig. 2 of [5] are mixed and inconclusive. 

The data on cc production coming from different experiments are not mutually quite com- 
patible but all fall into the broad band of theoretical prediction. This is possible because this 
band reflecting the current uncertainties of theoretical calculations is so broad that it can ac- 
commodate results differing by a factor of 2.5! The analogous band narrows somewhat for the 
bb production, but even there the lower and upper range of theoretical predictions differ by a 
factor of 2. In this case the two data points available are a factor of 3 above the median of 
theoretical predictions, but the experimental errors are too large to draw any definite conclu- 
sions. If confirmed, it would, however, be difficult to accommodate this result theoretically as 
one expects perturbative QCD to be better applicable to bb production than to the cc one. In 
the words of [3] "there is a serious discrepancy with bottom data". 

Theoretical "errors" come in part from the uncertainty in the values of m c and but also 
from the dependence of finite order QCD approximations on the renormalization and factor- 
ization scales and schemes. On the experimental side, major part of the error bars reflects the 
uncertainties in the extrapolation of visible cross section into the full phase space, which can be 
reduced only if better understanding of the theory is achieved. In such a situation it appears use- 
ful and timely to reanalyze the theoretical framework currently used for analyses of heavy quark 
production in 77 collisions, with particular attention to the renormalization and factorization 
scale invariance. This paper was motivated in part by nontrivial results of the detailed analysis 
of these dependencies in the case of heavy quark production in pp and pp collisions performed 
in [6]. Using this analysis as guidance, I will discuss the approximation currently used for the 
evaluation of (1(77 — > QQ) and point out its shortcomings. In particular, I will argue that it 
does not represent genuine next-to-leading order (NLO) QCD approximation and show that the 
missing ingredient are direct photon contributions of the order a 2 a 2 s . These terms, which have 
not yet been calculated, come in three classes, each playing its specific role. 

The paper is organized as follows. In the next two sections basic facts relevant for our 
discussion are recalled and the conventional treatment of (7(77 — > QQ) reviewed. In Section 
4 two distinct attitudes toward the meaning and content of the terms "leading" and "next-to- 
leading" orders of QCD are introduced and compared. Direct and resolved photon contributions 
to 17(77 — ¥ QQ) are analyzed in Sections 5 and 6 respectively, followed in Section 7 by the 
discussion of phenomenological consequences of the present analysis. 

2. Basic facts and formulae on the structure of the photon 

The factorization scale dependence of parton distribution functions (PDF) of the photon is 
determined by the system of coupled inhomogeneous evolution equations for quark singlet and 
nonsinglet and gluon distribution functions 
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The leading order splitting functions k q °\x) = (x 2 + (1 — x) 2 ) and P^'{x) are unique, while 
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all higher order ones k q ,kg , P}^ , j > 1 depend on the choice of the factorization scheme 
(FS). The equations (2.1-2.3) can be recast into evolution equations for qi(x,M),q i (x,M) and 



(0), 



(2.4) 
(2.5) 
(2.6) 



G(x,M) with inhomogeneous splitting functions kq^ = 3e 2 k q U) . The couplant a s depends on 
the renormalization scale fx and satisfies the renormalization group equation 
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where, in QCD with nj massless quark flavours, the first two coefficients, (3q = 11 — 2n//3 and 
f3\ = 102 — 38nj/3, are unique, while all the higher order ones are ambiguous. These non- 
unique coefficients, together with the boundary condition on the solution of (2.7), define the so 
called renormalization scheme (RS). The boundary condition is conveniently specified by the 
scale parameter Ars, which depends on the RS and corresponds to such a value of fi for which 
a s (fi) = 00. The couplant a s is thus actually a function of the ratio ju/Ars, which behaves for 
large /a as a s (///A) « 47r//?o ln(/_t 2 /A 2 ). For the sake of brevity I shall, however, drop the explicit 
specification of the dependence on the RS and write a s (fji) only. 

Provided perturbative expressions are summed to all orders, the results for physical quan- 
tities are independent of both the renormalization scale and scheme, but for any finite order 
approximation the numerical results do depend on \x and RS and the consistency of the theory 
merely requires that their variations with /x and RS are of higher order of a s than those taken 
into account. However, as the variation of a s (p) with \x (as well as with the RS) is proportional 
to a 2 , we must include at least first two consecutive nonzero powers of a s in perturbative ex- 
pansions of physical quantities for the necessary cancellation to operate and for the finite order 
approximation to be performed in a well-defined RS of the couplant a s . 

Due to the presence of the inhomogeneous splitting functions k q and kc a general solution 
of the evolution equations (2.1-2.3) can be split into the particular solutions of the full inho- 
mogeneous equations and a general solutions, called hadron-like (HAD), of the corresponding 
homogeneous ones. A subset of the former resulting from the resummation of contributions 
of diagrams in Fig. 1, describing multiple parton emissions off the primary pure QED vertex 
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Figure 1: Diagrams defining the point-like parts of nonsinglet quark and gluon distribution functions. 
The resummation involves integration over parton virtualities r < M 2 and is represented by the junction 
of the blob and the 7 — > qq vertex. Partons going into the hard collision are denoted by x, t. 

7 — > qq and vanishing at M = Mq, are called point-like (PL) solutions. Due to the arbitrariness 
in the choice of Mq the separation of PDF into their point-like and hadron-like parts is, however, 
ambiguous. In general we can thus write (D = q, q, G) [7] 

D(x, M) = D PL (x, M, M ) + D RAD (x, M, M ). (2.8) 

The hadron-like parts can be represented by the same solid blobs as those of hadrons. Appro- 
priate graphical representation of the point-like ones, proposed in [8, 9], is shown in Fig. 1. By 
joining the primary 7 — > qq vertex to the solid blob it reflects its perturbative origin. 

The different nature of the UV renormalization of QCD quantities, generating the renormal- 
ization scale and scheme dependence of a s , and IR "renormalization" involved in the definition 
of PDF represents the main reason for keeping the factorization and renormalization scales M 
and ji as independent free parameters. The former sets the upper bound on the virtualities of 
partons included in the definition of PDF and comes from IR part of loop corrections as well 
as integrals over real parton emissions. The renormalization scale, on the other hand, defines 
the lower bound on virtualities included in the renormalized colour charges, masses and fields 
and comes entirely from loops. There is no compelling theoretical reason why these two scales 
should be identified [10]. 

3. (7(77 — > QQ): the conventional approach 

In [1, 2, 3] the "next-to-leading order" QCD approximation is defined by taking the first two 
terms in expansions of direct, as well as single and double resolved photon contributions 

o- diT = <r§> + <r$>a s (jji) + 4J(M, fi)a 2 s (fi) + <r®(M, ^ s (jjl) + ■■■, (3.1) 
a ST = a^\M)a s ( f i) + a^(M, f i)a 2 s ( j i) + a^(M, f i)a 3 s ( f i) + ---, (3.2) 

a dr = 4 2 r ) (M)a2( /u)+fT ( 3 r )(M,/,)a3(^) + ... (3.3) 

to the total cross section 17(77 — * QQ) f° r the heavy quark pair production in 77 collisions 

0-(77 -»• QQ) = O-dir + + C7 dr . (3.4) 

In [1, 2, 3] three light quarks (u, d, s) were considered as intrinsic for the evaluation of (1(77 — > cc) 
and four (u, d, s, c) in the case of bb production, which is appropriate taking into account that 



3 



mi, S> m c . Although the general strategy of the GRV group is to consider as intrinsic quarks (of 
both the photon and proton) the light quarks u, d, s only, the GRV set used in [1] parameterizes 
the effects of heavy quark production far above their thresholds by means of the charm and 
bottom quark distribution functions and is thus compatible with [1]. 

Note that while af^ and are functions of s/ttiq only, the lowest order single and double 

resolved photon contributions and cr^) depend, via PDF of the light quarks and gluons, also 
on the factorization scale M, assumed for simplicity to be the same for both colliding photons. 
All higher order coefficients in (3.1-3.3) depend on both the factorization and renormalization 
scales M and ji. The lowest order term in (3.1) comes from pure QED and equals 



For the direct photon contribution (3.1) the renormalization and factorization scales \x and M 

f<y\ 

appear first in <r^, which, however, is not included in [1, 2, 3]. 

Defined in this way, the direct, single resolved and double resolved contributions start and 
end in (3.1-3.3) at different powers of a s . In the [1, 2, 3] this is justified by claiming that PDF of 
the photon behave as a/a s , and consequently all three expansions (3.1-3.3) start and end at the 
same powers (a s )° = 1 and a s , respectively. However, as argued in detail in [11], the logarithm 
InM 2 characterizing the large M behaviour of PDF of the photon cannot be interpreted as 
l/a s (M) , as it comes from integration over the transverse degree of freedom of the purely QED 
vertex 7 — ► qq\ If QCD is switched off by sending, for fixed Mq, A — > 0, quark and gluon 
distribution functions of the photon approach their finite QED expressions 



qi (x,M) -+ q? EB (x,M) = -^Seffef (x)ln-2, G{x,M) -+ G^ ED (x,M) = 0, (3.6) 



where Mq regularizes the parallel singularity coming from the 7 — > qq vertex by setting the lower 
limit on the virtuality of the quark or antiquark going to the hard collision. The relation (3.6) 
then implies that the single and double resolved photon contributions (3.2,3.3) vanish in this 
limit, because a s (fi) does so. In the absence of QCD we thus get, as we must, the purely QED 
result cr^. Had the PDF of the photon really behaved as a/a s , we would, on the other hand, 
expect finite contributions from the lowest order single and double resolved photon contributions 
(3.2,3.3) even in the limit of switching QCD off. I will discuss this limit, together with the large 
M behaviour of PDF of the photon in detail in subsection 6.1. Taking into account the separation 
(2.8) we can distinguish 5 classes of resolved photon contributions: 

• Single resolved photon using 

- hadron-like parts of PDF (cr sr h), 

- point-like parts of PDF (cr srp ). 

• Double resolved photon using 

- hadron-like parts of PDF on both sides (<7drhh)> 




(3.5) 



where s denotes the square of total 77 collision energy, f3 = 
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— hadron-like parts of PDF on one side and point-like ones on the other (<7drhp)> 

- point-like parts of PDF on both sides (odrpp)- 

With this subdivision of a SI and <7d r in mind we can rewrite (3.4) as follows 

(7(77 — > QQ) = (Tdir + Csrh + Csrp + Cdrhh + CTdrhp + Cdrpp- (3.7) 

4. Discourse on semantics: denning LO and NLO in QCD 

Although the definition of the concepts "leading-order" (LO) and "next-to-leading-order" (NLO) 
QCD approximations is primarily a matter of semantics, it is in my view preferable to use the 
terminology that respects the basic fact that we have in mind orders of QCD, rather than the 
total number of terms taken into account in expansions (3.1-3.3). Counting only the powers 
of a s will allow us to associate the term "NLO" with calculations performed in a well-defined 
renormalization scheme of the couplant ct s as well as factorization scheme of PDF. In this context 
it is useful to recall the meaning of theses terms for three physical quantities related to (3.4). 

4.1 <r(e + e - — > hadrons) 

For proper treatment of the direct photon contribution a^ T in (3.1), the total cross section for 
e + e~ annihilations into hadrons at the cms energy Q provides particularly suitable guidance. 
For n f effectively massless quark flavors we have 

CThad(Q) = 4ad(Q) + a,Gu)4ad(Q) + «s (rfad (QM + " " " = ^ad^X 1 + KQ)), (4-1) 

where the lowest order term o~^ d (Q) = (4-ira 2 / Q 2 ) e 2 comes, similarly as a^j. in (3.1), from 
pure QED, namely the diagram in Fig. 2a, whereas genuine QCD effects are contained in 



r(Q) = °>M 



TT 



1 H ri(Q/n) + • 
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(4.2) 



The lowest order contribution to r(Q) results from the sum of the integral over the cross section 
for real gluon emissions (Fig. 2b) with the interference term between the basic QED diagram of 
Fig. 2a and one loop corrections to it, shown in Fig. 2c. The next-to-lowest order contribution to 
(4.2), (a s /n) 2 ri, results from summing the integral over the double parton emissions, exemplified 
by the diagrams in Fig. 2d, with the integral over the interference term between single gluon 
emission (Fig. 2b) and one loop corrections to it (for examples see Fig. 2e) , and the interference 
term between the QED contribution of Fig. 2a and the two loop corrections to it (not shown). 
At each order of a s the finiteness of the sum is guaranteed by the KLN theorem. 

For the quantity (4.1) nobody calls the lowest order term cr[^ d (Q) the "leading-order" and 
the next one, i.e. o-^ d (Q)a s ((i) / tt , the "next-to-leading order" QCD approximations, but these 
terms are applied to genuine QCD effects contained in r(Q). To work in a well-defined renormal- 
ization scheme of a s requires including in (4.2) at least first two consecutive nonzero powers of 
a s (n), because only the coefficient r\(Q / fi) /it 2 standing by a 2 depends on the renormalization 
scale and scheme. The renormalization scale invariance of the NLO approximation to tThad(Q) 
implies that the implicit ^-dependence of the leading order term a s (/x)/7r in (4.2) is cancelled to 
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Figure 2: Feynman diagrams contributing to the cross section a ee (Q) up to the order a 2 a 2 8 . Diagrams 
with real gluon emissions off and loop corrections on the upper quark leg are understood to be included. 

order a 2 by explicit //-dependence of the coefficient r\(Q/ /j,)/it 2 . The In fj, 2 term necessary for 
this cancellation comes from loop corrections described by diagrams like those in Fig. 2e. Note 
that the loops in Fig. 2c contribute to the renormalization of a, rather than of a s . Because 
a s (fi) is a monotonous function of \i spanning the whole interval (0,oo), the inclusion of first 
two consecutive nonzero powers of a s is also a prerequisite for the applicability of any of the 
scale fixing methods [12, 13, 14] available on the market. Without it, there is no preferred scale 
and, consequently, the LO QCD approximations become entirely arbitrary. 

For purely perturbative quantities the association of the term "NLO QCD approximation" 
with a well-defined renormalization scheme is a generally accepted convention, in my view worth 
retaining for physical quantities in any hard scattering process. 

4.2 F$(x,Q 2 ) 

For physical quantities involving PDF proper definition of the terms "leading" and "next-to- 
leading" order concerns beside the renormalization scale and scheme of a s also the factorization 
scale and scheme of these PDF. 

For proton the structure functions F^x^Q 2 ) is given as the convolution of PDF with the 
corresponding coefficient functions. Dropping for simplicity quark charges we generically have 

F$(x, Q 2 ) = q(M) ® C q (Q/M) + G(M) ® G G (Q/M), (4.3) 

where quark and gluon distribution functions satisfy the homogeneous evolution equations and 

C q (x,Q/M) = 5(l-x) + ^C( 1 \x,Q/M) + (^YcP(x,Q/M,Q/ri + .-. (4.4) 

C G (x, Q/M) = + ^C£\x, Q/M) + 0^-) 2 Cg\x, Q/M, Q/fi) + • • • . (4.5) 

The standard definition of LO QCD approximation to F 2 p (x, Q 2 ) amounts to setting 

Fl LO (x,Q 2 )=q(x,M), (4.6) 

where M oc Q and q(x,M) satisfy (2.1- 2.3) with P^p on their r.h.s. only. At the NLO 
the terms proportional to P^\Cq^ and Cq* are included as well. All these functions depend 
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on the factorization scheme, and the latter two also on M, and must therefore be included 
in any analysis which aims at working in a well-defined such scheme. Although i ? |' L Q(x,Q 2 ) 
involves only the first, purely QED term in (4.4), it describes QCD effects, because they drive 
the M-dependence of q(x,M). At the NLO, the factorization scale dependence of q(x,M) 
in (4.3) is cancelled by the explicit M-dependence of the coefficient C^\Q/M), similarly as 
renormalization scale dependence of a s (n) in the first term of (4.2) is cancelled to the order 
a 2 by the jU-dependence of the coefficient ri(Q/fi). The standard definitions of LO and NLO 
approximations for F^ix^Q 2 ) thus conform to the basic requirement that the term "leading 
order" describes the lowest order contribution affected by QCD effects and that working in a 
well-defined renormalization as well as factorization schemes requires at least the next-to-leading 
order approximation. 

4.3 F^x,Q 2 ) 

The conventional analysis of photon structure function employs the same definition of the terms 
"leading" and "next-to-leading" orders as that used in [1, 2, 3]. I have discussed its shortcomings 
in [11] and will therefore merely recall the main points relevant for this paper. For the real photon 
F^(x, Q 2 ) is given as the sum of convolutions 

-F](x, Q 2 ) = q m (M) ® C q {Q/M) + ^-<5nsC 7 + 

X ZTT 

(e 2 )£(M) ® C q {Q/M) + ^( e 2 )5 s C 7 + (e 2 )G(M) ® C G (Q/M) (4.7) 

Z7T 



of PDF of the photon with the coefficient functions C q and Cq, given in (4.4-4.5), and 

C 7 (x, Q/M) = Cf(x, Q/M) + S!M C W( X , Q/M) + (^jjp) C^{x, Q/M) + • • • , (4.8) 
where the lowest order term in (4.8) 

C(°\x,Q/M) = (x 2 + {l-x) 2 ) 



, M z , l-x 



+ 8x(l-x)-l (4.9) 



comes, similarly as k q in (2.4) and cr^ d in (4.1), from pure QED. Throughout this subsection 
I will discuss only the point-like part of F^XjQ 2 ), which is related to the point-like parts of 
PDF of the photon. The QED contribution to F^(x,Q 2 ) 

-, n f 

^ ED (x,Q 2 ) = £e 2 (q? ED (x,M) +^ ED (x,M)) + ^-6n f (e*)CW(x,Q/M) (4.10) 
i=i 

comes from the diagram in Fig. 3a and is actually independent of M because the M-dependence 
of q® ED (x,M) is cancelled by that of C^°\x,Q/M). Switching on QCD implies that terms 
proportional to nonzero powers of a s in the splitting as well as coefficient functions (2.4-2.6) 
and (4.4,4.5,4.8) are taken into account. 

It appears natural to define as the leading order QCD approximation to F^ the one that 
takes into account those of the above terms proportional to a s which contribute to the 0(a s ) 
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a) 
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Figure 3: Examples of diagrams contributing to F% in QED (a), 0(a s ) (b-c) and 0{a 2 s ) (d-i) QCD. 



correction to (4.10). These are k^ and P^> in (2.4-2.6) and C q y> and C^' in (4.4) and (4.8). 
The remaining 0(a s ) terms, k^ and Cq\ do not contribute to a s correction to because 
G^ ED (x,M) = 0. The functions kq~\ as well as result from evaluation of the order 
a 2 a s diagrams 1 , exemplified by those in Fig. 3b, c. It is worth emphasizing that despite the 
presence of loops involving gluons, no renormalization of a s is performed at this order, because 
the loops, like that in Fig. 3b, contribute to the UV renormalization of a, not a s l In this 
definition the LO QCD approximation to thus involves (see [11] for explicit formula) in 
addition to the purely QED quantities k g °^ and also kq~\ and Note that the 

conventional LO approximation takes into account out of these 6 quantities only two of them: 
4 0) and P^. 

According to the approach advocated in [11], the next-to- leading order QCD approximation 
to F2 includes in addition to the quantities mentioned above also the splitting and coefficient 
functions proportional to a 2 , i.e. k q 2 \ P[j \ Cq 2 \ and takes into account also the contribution 

generated by the gluonic splitting and coefficient functions k^ and Cq\ The evaluation of 
these quantities involves order a 2 a 2 s diagrams, exemplified by those in Fig. 3d-i. These include 
loop corrections to the LO diagrams (like those in Fig. 3e-g) as well as two parton emissions, 
exemplified by diagrams in Fig. 3d,h-i. At this order the renormalization of a s starts to operate 
by removing the UV divergencies coming from loops in Fig. 3e-g (and others). Note that the 
conventional NLO approximation to F^ (see, for instance, [15]) includes k q °\ k^\ Pjj \ Pij \ ^7°^ 
and Cq 1 ^ , but not k q 2 ^ , C q 2 ^ , , . This is particularly intriguing in the case of the coefficient 



>(o) ■ 



(i) 



function which stands in (4.8) by power aa s and is actually known! As k q *' ,C q *' ,Cq' and 
Cj 2 \ have not yet been calculated, a complete NLO QCD approximation of F^ defined in the 
sense advocated in [11], cannot at the moment be constructed. 

The fact that contrary to the case of F^, C q ^ enters F^ first at the NLO is due to the fact 
that there is no point-like, perturbatively calculable coupling of the proton to quarks (or gluons) 
that would generate the inhomogeneous splitting function analogous to k q °^ , and thus the pure 
QED contribution analogous to (3.6). 



(2) n {2) n (l) 



x By "order of diagram" I mean in the case of tree diagrams the order of its square and in the case of diagrams 
with a loop the order of the interference term of this diagram with the corresponding tree one(s). 
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5. Direct photon contribution to (7(77 — > QQ) 

In [1, 2, 3] the purely QED contribution is considered as the leading-order and the sum 

as the NLO approximation of the direct photon contribution Odir- Note that (5.1) is of the 
same form as the sum of first two terms in (4.1) and a s (p) is thus the only place where \x 
appears. Consequently, cannot be associated to a well-defined renormalization scheme of 
a s and therefore does not deserve the label "NLO" even if the NLO expression for a s (n) is used 
therein. For QCD analysis of o^ir m a well-defined renormalization scheme the incorporation 
of the third term in (3.1), proportional to a 2 a 2 , is indispensable. We need in particular the 
diagrams like those in Fig. 4j-k, which involve loops contributing to the renormalization of a s . 
The regular, /i-dependent part of their contributions provides the In fi 2 term canceling the \i 
dependence of a s (jx) in the second term of (5.1). 




Figure 4: Examples of diagrams contributing to (Tdi r in QED (a), leading (b-d) and next-to-leading (e-1) 
orders of QCD. Heavy quark are denoted by dashed, light quarks by solid lines. 

At the order a 2 a 2 diagrams with light quarks appear as well and we therefore distinguish 
three classes of contributions, differing by the charge factor CF (e<g and e q denoting charges of 
heavy and light quarks respectively): 

Class A: CF = eg. Comes from diagrams in which both primary photons couple to heavy 
quarks or antiquarks. The diagrams of this class may, as those in Fig. 4i-j, contain also 
light quarks. As emphasized above the diagrams containing loops, like those in Fig. 4j- 
k, are vital for the renormalization of a s . Despite the presence of mass singularities of 
individual diagrams coming from gluons and light quarks in the final state and loops, the 
KLN theorem guarantees that at each order of a s the sum of all contributions of this class 
to o"dir is finite. 

Class B: CF = ege 2 . Comes from diagrams in which one of the primary photons couples to a 
heavy and the other to a light quark-antiquark pair, like that in Fig. 4h. For massless light 
quarks this diagram has initial state mass singularity, coming from the region of vanishing 
light quark virtuality T q , which is removed by introducing the concept of the resolved 
photon. This implies subtracting from the corresponding cross section the integral over 
the double pole l/(T g Tc) and putting it into the point-like part of the gluon distribution 
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function appearing in the lowest order single resolved photon diagram of Fig. 5a. Similarly, 
the single pole term l/r g is subtracted and put into the light quark distribution function 
entering the next-to-lowest order single resolved photon diagram if Fig. 5c. 

Class C: CF = ei. Comes from diagrams in which both photons couple to light qq pairs, as 
those in Fig. 41. In this case the analogous subtraction procedure relates it to the lowest 
order double resolved photon contribution of the diagram in Fig. 5h. 

Because of different charge factors, the classes A, B and C do not mix under renormalization 
of a s and factorization of mass singularities. Such mixing does, however, occur, within each 
of these classes. At the order a 2 a 2 s all three classes of contributions are needed for theoretical 
consistency, albeit each for different reason. The class A is needed if the calculation is to be 
performed in a well-defined RS. In the next Section we shall see that classes B and C must 
be taken into account for factorization scale invariance of the sum (3.4) of direct and resolved 
photon contributions. 



6. Resolved photon contribution to (7(77 — > QQ) 

Because the factorization mechanism plays crucial role in my arguments let me recall how it 
works for heavy quark production in pp collisions, where the NLO approximation for a(pp — > 
QQ) involves convolutions of PDF with partonic cross sections up to the order <x? . Schematically 



cr 1N|ijU (pp — > QQ) = D\(M) <g) [a 2 s (fi)a^ + a 3 s (n)a^ (M, fi)^ ® D 2 (M), (6.1) 

where PDF D\(M) and Di^M) of the colliding protons satisfy the homogeneous evolution 
equations with first two terms in (2.6) taken into account. Factorization scale invariance of (6.1) 
is guaranteed by the fact that the M-dependence of D\(M) and Z?2(M) in the lowest order 
term of (6.1) is cancelled to the order considered by the M-dependence of (M, //). Graphical 
representation of this cancellation mechanism exploits the fact that the homogeneous part of the 
evolution equations (2.1-2.3) relate by what I call homogeneous factorization a given diagram 
with partonic cross section at the order a 2 , with two diagrams (with incoming quark and gluon 
respectively) at the order c? s and higher. For the NLO approximation (6.1) only the lowest order 
splitting functions P^p appear in this cancellation. 

In [6] detailed analysis of the renormalization and factorization scale dependence of (6.1) 
has been performed under the assumption fi = M. The results, summarized in Figs. 13-15 of 
[6], show that the sensitivity of cr NLO (pp — ► QQ) to the variation of M depends sensitively on 
the ratio ^/s/rng, as well as on PDF of colliding particles. They also demonstrate that contrary 
to the traditional expectation the NLO QCD approximation is not necessarily less dependent 
on the factorization scale M than the LO ones ! Whereas at y/l = 62 GeV a NLO (pp -► bb) 
is a weaker function of M than d LO (pp — > bb), and has even a stationary point, this is not 
true at yfs = 630 GeV. There the NLO prediction is a monotonous function of M, which is 
even steeper than the LO one! Similar lesson follows from comparing different processes at 
the same energy: whereas <7 NLO (pp — > bb) has, as mentioned, a stationary point for ^fs = 62 
GeV, cx N (pp — ► bb) is again a monotonous function of M. To assess the stability of QCD 
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calculations of 17(77 ~^ QQ)i quantitative analysis of their factorization and renormalization 
scale dependence along the lines of [6] would be very helpful, better still with fj, and M treated 
as separate free parameters. 



6.1 Single resolved photon: the point-like part 

For the point-like part of PDF of the photon the presence of the inhomogeneous splitting function 
in the evolution equations (2.1,2.3) implies additional relation, which I call inhomogeneous 
factorization, to distinguish it from the homogeneous one introduced above, which corresponds 
to P^p ■ For instance, the lowest order single resolved photon diagram of Fig. 5a is related 
by homogeneous factorization to the single resolved photon diagram of Fig. 5c. Both of these 
diagrams are included in the approximation 



o-W(M,v) ee a^JM)a s ^) + 4 2 T l(M,fi)a 2 M 



(6.2) 



employed in [1, 2, 3]. The latter diagram is related via the inhomogeneous factorization to the 
direct photon diagram of Fig. 5b, which is of the order a 2 a 2 s and thus is not included in the 
approximation (5.1) of er dir . 



\ . \ . v.. v.... v.... \ 





Figure 5: Some of Feynman diagrams contributing to er srp in the leading (a) and next-to-leading (c-f) 
orders of QCD, together with the related direct photon diagrams of the order o?o? s (b,g) and lowest order 
double resolved photon diagram (h). 

The inhomogeneous factorization, which connects direct and resolved photon diagrams at 
the same order of a s , has a simple interpretation, reflecting its role in removing the mass 
singularities of direct photon contributions associated with the point-like coupling of initial 
photons to light quarks. For the direct photon diagram in Fig. 5b (which coincides with the 

diagram in Fig. 4f) this procedure has been outlined already in Section 5. Let me recall that its 

(2) 

finite, M-dependent part contributes to class B part of a d ^ and cancels the dominant part of 
the factorization scale dependence of single resolved photon diagram in Fig. 5c. The inclusion 
of class B direct photon contributions of the order a 2 a 2 is thus vital for factorization scale 
invariance of the sum (3.4) of direct and resolved photon contributions. In similar way, the 
factorization scale invariance implies that single resolved photon diagram of Fig. 5f must be 
considered together with class C direct photon diagram of the order a 2 a 2 , shown in Fig. 5g. 
This diagram is in turn related by inhomogeneous factorization to the point-like part of lowest 
order double resolved photon diagram of Fig. 5h, which is of the same order a 2 a 2 and will be 
discussed below. These considerations apparently contradict the statement in [1] 
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In contrast to heavy-quark production in hadron collisions we expect small variation 
of the resolved-^/ cross section with [i at the Born level, since the \i fall off' of the 
parton cross section [~ a k {n)] is neutralized asymptotically by the increasing number 
of gluon and quark partons in the photons ~ [aj 1 (fi)] k , k = 1 or k = 2 for 1- or 
2-resolved 7. 

I will now argue that the asymptotic behaviour of Born terms in the resolved photon contribu- 
tions, mentioned in the above quotation, does not actually imply factorization scale invariance 
of the approximation used in [1, 2, 3]. The large M behaviour of the solutions of the evolution 
equation (2.1-2.3) is the same for all three PDF and reads (in momentum space) 



q NS (n,M) = _5 NS fcf(n) 



£(n,M) = — (n) 



a 



G(n,M) = -^f(n) 



InM 2 , 



1 - P GG {n) 



L(l - P qq (n))(l - P GG (n)) - P qG (n)P Gq (n) 

P G g(n) 

;i " P gq (n))(l ~ P GG (n)) - P qG (n)P Gq {n) 



E(n, M) 



PqM 

1 - P GG (n) 



(6.3) 

InM 2 , (6.4) 
InM 2 , (6.5) 



where Pij(n) = 2P^ (n) /(3q . Although gluons must be radiated off the primary quarks, which 
costs powers of a s , the InM 2 rise of point-like part of quark distribution function drives the 
logarithmic rise of G(n, M) as well. Consequently, aill(M)a s (M), the lowest order term in (6.2) 
(where we set /U = M) approaches at large M a constant, whereas in hadronic collisions the first 
term in (6.1) vanishes as M — > 00. This difference might seem important, but actually is not, 
because the factorization scale invariance concerns the variation of finite order approximations 
with M, not their magnitude! Moreover, the values of M used in [1], i.e. M = \/2mQ is far 
from asymptotic in any case. 

The relations (6.3-6.5) determine the asymptotic behaviour of PDF of the photon for fixed 
value of the QCD parameter A. The claim that PDF of the photon behave as a/a s {M) must 
therefore be interpreted merely as a shorthand for this large M behaviour at fixed A. Interpreting 
InM 2 as InM 2 /A 2 oc a s (M), would lead us to obviously wrong conclusion that PDF of the 
photon blow up to infinity when we switch QCD off by sending A — ► 0! On the other hand, as 
(6.3-6.5) do not explicitly contain A, they seem to hold for any value of A and, consequently, 
one might be tempted to conclude that they survive the limit A — > as well. That, however, is 
not the case. 

To find what happens with the asymptotics (6.3-6.5) when we switch QCD off we must 
reverse the order of limits and first send A — ► for fixed finite Mo and M. Doing this we get, 
of course, the pure QED formulae (3.6), which differ dramatically for the gluon and quantita- 
tively (i.e. by the presence of the square brackets in (6.3-6.4) also for the quark distribution 
functions. The corresponding asymptotic behavior of quark and gluon distribution functions, 
can be obtained formally also directly from (6.3-6.5) by sending (3q — ► 00, which implies setting 
Pij = 0, Vi,j there. This reflects the fact that the factor 2//3q multiplying in (6.3-6.5) all the 
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splitting functions results from the following limit 

a s (M) 2 2 InM 2 2 , N 

lim -l^lnM 2 = — lim - = 6.6 

M^oo 2tt (3 M^oo ln(M 2 /A 2 ) p 

Sending a s — ► for fixed M can be achieved either by sending A — > for fixed finite M, or, 

formally, for any M by sending (3q — ► oo. The latter option is applicable even for the asymptotics 

(6.3-6.5). Clearly, the order of the limits A —> and M —> oo matters. 

There is also a second difference between QQ production in pp and 77 collisions. At the 

NLO the q PL contributes (see diagram in Fig. 5c) to the second term in (6.2), as it does the 

quark distribution functions of the proton to the second term in (6.1). But whereas in the latter 

case the variation of this term with M is of one order of a s higher than that taken into account 

in (6.1), for 77 collisions the presence of kf^ implies that the variation of a^la 2 s {M) with M is of 

the same order as this contribution itself. As argued above, the cancellation of this dependence 

requires the inclusion of direct photon contributions of class B and order a 2 a 2 . 

(12) 

In summary, the approximation <7s rp (M, /j) defined in (6.2) has the properties of genuine 
NLO QCD approximation with respect to the renormalization of a s , but without the inclusion 
of direct photon contributions of classes B and C, like those in Fig. 5b, g, it is not factorization 
scale invariant to the order considered. 

6.2 Single resolved photon: the hadron-like part 

This part of (1(77 — ► QQ) has the same structure as the direct photon contribution to the cross 
section for heavy quark production in 7p collisions and I will therefore discuss the main features 
of the approximation 

^h 2) ( M ^) = ^l(M)a s (») + a£l(M,»)a 2 M. (6.7) 

(23) 

in the next subsection together with those of o"^ rhp , which plays the role of the resolved photon 
contribution in 7p collisions. 

6.3 Double resolved photon: the hadron-like point-like part 

This part of double resolved photon contribution used in [1, 2, 3] 

A = 4 2 rh P ( M )«^) + 4rhp( M > M)c&M) (6-8) 

has the same properties as the point-like part of resolved photon contribution to cr(7P — > QQ). 
In fact, only the first term in (6.8) needs to be included in the sum 

^'W^+CplM.p) (6.9) 

to make it complete NLO QCD approximation as far as both the renormalization of a s and 

(2) 

the definition of PDF of the photon are concerned. Note that the M-dependence of cr^^M, /j,) 

induced by the inhomogeneous splitting function compensates the M-dependence of the part 
(2) 

of *£i(M,fj) that remains after the subtraction of the singular term coming from the vertex 
7 -> qq. 
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6.4 Double resolved photon: the hadron-like hadron-like part 

This part of double resolved photon contribution has the same structure and properties as (6.1) 
and, consequently, the sum of first two terms in (3.3) 

*!£L(M, fj) = cr^ hh (M)a 2 M + 4 3 2 h (M, »)a 3 M (6.10) 

has all the features of genuine NLO QCD approximation, both as far as the renormalization of 
a s and the definition of PDF are concerned. 

6.5 Double resolved photon: the point-like point-like part 

The lowest order contribution to this part of (1(77 — > QQ) comes from the convolution of 
point-like parts of PDF of both photons with a 2 cross sections of qq — > QQ (see Fig. 5h) and 
GG — > QQ (not shown) partonic subprocesses. The former contribution is needed to render the 
point-like part of single resolved photon contribution (6.2) coming from the diagram in Fig. 5f 
factorization scale invariant. But as argued in subsection 6.1, the diagram in Fig. 5h alone is 
not sufficient for this purpose and class C direct photon diagram of Fig. 5g must be included 
as well. For the double resolved photon contribution the same diagram is needed for theoretical 
consistency already at the lowest order, i.e. for the cancellation of part of the M-dependence of 

7. Phenomenological consequences 

In the preceding sections I have shown that the approximation used in [1, 2, 3] 

a DKF (M, n) = ( M ) + c4 12 ) (M, M ) + a£ 3) (M, M ), (7.1) 

based on the sum of first two term in each of the expansions (3.1-3.3), does not represent genuine 
NLO QCD approximation despite the fact that it contains partonic cross sections up to order 
and uses the NLO form of a s {fj). The missing direct photon contributions of the order a 2 a 2 s 
come in three classes, which differ by the overall charge factor, reflecting the presence or absence 
of light quark pairs and the way they couple to the primary photons. Each of these three classes 
is needed for different reason, but all are needed to make the expression 

cx NLO (M, fj,) = af 2) (M, fx) + (M, fi) + (M, M ) (7.2) 

genuine NLO QCD approximation. Without them, (7.2) reduces to (7.1), which is complete and 
theoretically consistent merely to the leading order of QCD, both as far as the renormalization 
of a s and the factorization of mass singularities into PDF are concerned. 

As higher order calculations involving heavy quarks are difficult to perform, and complete 

(2) 

calculation of a A ^ thus not in sight, it is important to find ways of estimating their numerical 
importance. 
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7.1 Class A 

This class of direct photon contributions, related to the renormalization scale and scheme depen- 
dence of (Jdir, does not mix with the resolved photon contributions at any order and can therefore 
be considered separately. Because the approximation (5.1) is a linear function of a s (fx), one can 
by choosing appropriate value of jJL get arbitrary value of a s {pL) and thus also of (5.1), even if the 
NLO form of the //-dependence of a s ([i) is used. The existence of a well-defined "natural" scale, 
like tuq in our case, is of no help in this respect, as a s (n) depends beside the renormalization 
scale fi also on the renormalization scheme RS, and at the LO there is no criteria for selecting 
the "natural" RS. Without the inclusion of the direct photon contributions of this class and 
order a 2 a 2 , there is no way how to estimate the importance of higher corrections. According to 
[1], cF^ T oi s (fi) gives, setting = M = \^2mQ, about 25 — 35% of the pure QED contribution <r^. 
This is a sizable correction, which can easily be doubled by inclusion of higher order corrections 
of this class. 

7.2 Classes B and C: controlling the scale ambiguities 

These classes are required to render the sum of point-like parts of single and double resolved 
photon contributions factorization scale invariant. As none of them is included in [1, 2, 3], the 
sum of resolved photon contributions 

a^ F (M, fi) ee (M, n) + af ] (M, fi) (7.3) 

used in [1, 2, 3] is expected to exhibit monotonous dependence on the factorization scale M. 
On the other hand, as the point-like parts of PDF of the photon can be considered as a way of 
approximate evaluation of higher order perturbative corrections [7], one may wonder whether 
by clever choice of the factorization scale M one could approximately include direct photon 
contributions of the classes B and C in the point-like parts of single and double resolved photon 
ones. To establish whether this can be arranged and how, requires, however, a quantitative 
analysis of factorization and renormalization scale dependence of the expression (7.3). Specif- 
ically, one should plot for given s 77 and mg, the quantity (7.3) as a two-dimensional function 
of ji and M. Such a plot would give us a clear idea of the stability of (7.3) with respect to the 
variations of [i and M and, possibly, suggest better choices of these scales than those adopted 
in [1] (/x = M = \f2m Q ) or [2, 3] (M = 2/i = 2m T ). 

8. Summary and conclusions 

The fact that the conventional NLO analyses of heavy quark production in 77 collisions do 
not include direct photon contributions of the order a 2 a 2 s represents a serious shortcoming 
preventing us from drawing definite conclusions from the comparison of existing data on cc and 
bb production with currently available QCD calculations. 

The missing direct photon contributions of the order a 2 a 2 come in three classes, depending 
on the overall charge factor. One of them, proportional to eg, concerns exclusively the direct 
photon contribution to (7(77 — > QQ), and is vital for establishing the genuine NLO character of 
o"dir- The other two classes of direct photon contributions of the order a 2 a 2 are needed for the 
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factorization scale invariance of the sum of direct and resolved photon contributions. To assess 
their numerical relevance, quantitative investigation of the factorization and renormalization 
scale dependence of the approximation (7.3) would be very helpful. 

I am grateful to M. Kramer for correspondence concerning the treatment of intrinsic charm 
in the bb production. 

References 

[1] M. Drees, M. Kramer, J. Zunft, P. Zcrwas, Phys. Lett. B 306 (1993) 306. 

[2] M. Kramer, E. Lacncn, Phys. Lett. B 371 (1996) 303. 

[3] S. Frixionc, M. Kramer, E. Laenen, Nucl. Phys. B 571 (2000) 169. 

[4] D. Buskulic ct al. (ALEPH Collab.), Phys. Lett. B 355 (1995) 595. 
D. Buskulic et al. (ALEPH Collab.), ALEPH 2000-70 
G. Abbiendi et al (OPAL Collab.), Eur. Phys. J. C16 (2000) 579 
C. Abbiendi ct al. (OPAL Collab.), OPAL Physics Note 455 

A. Csilling, Proceedings of PHOTON2000, Ambleside, England, September 2000, ed. A. Finch 
M. Acciarri et al. (L3 Collab.), Phys. Lett. B 453 (1999) 83 
M. Acciarri et al. (L3 Collab.), Phys. Lett. B 467 (2000) 137 
M. Acciarri et al. (L3 Collab.), Phys. Lett. B 503 (2001) 10 

[5] E. Laenen, S. Frixionc, M. Kramer, Proceedings of PHOTON2000, Ambleside, England, September 
2000, cd. A. Finch 

[6] G. Altarelli, M. Diemoz, G. Martinelli, P. Nason, Nucl. Phys. B 308 (1988) 724 

[7] J. Chyla, M. Tasevsky, Phys. Rev. D62, (2000) 114025 

[8] C. Friberg, T. Sjostrand, JHEP09(2000)010 

[9] J. Chyla, hep-ph/0102100 

[10] P. Aurenchc, R. Baier, M. Fontannaz, D. Schiff, Nucl. Phys. B 286 (1996) 509 

[11] J. Chyla, JHEP04(2000)007 

[12] P. M. Stevenson, Phys. Rev. D 23 (1981) 2916 

[13] G. Grunberg, Phys. Rev. D 29 (1984) 2315 

[14] S. Brodsky, G. Lepage, P. Mackenzie Phys. Rev. D 28 (1983) 228. 

[15] M. Gliick, E. Reya, I. Schienbein, Phys. Rev. D60 (1999), 054019 



16 



